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This paper presents an analytical study of the coexistence of different transport regimes in quasi- 
one-dimensional surface-disordered waveguides (or electron conductors). To elucidate main features 
of surface scattering, the case of two open modes (channels) is considered in great detail. Main 
attention is paid to the transmission in dependence on various parameters of the model with two 
types of rough-surface profiles (symmetric and antisymmetric). It is shown that depending on the 
symmetry, basic mechanisms of scattering can be either enhanced or suppressed. As a consequence, 
different transport regimes can be realized. Specifically, in the waveguide with symmetric rough 
boundaries, there are ballistic, localized and coexistence transport regimes. In the waveguide with 
antisymmetric roughness of lateral walls, another regime of the diffusive transport can arise. Our 
study allows to reveal the role of the so-called square-gradient scattering which is typically neglected 
in literature, however, can give a strong impact to the transmission. 
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I. INTRODUCTION 

The transport of quantum quasi-particlcs through thin 
metal films and semiconductor nano-structures, such 
as nanowires and strips, supperlattices and quantum- 
well systems, has been a hot topic since more than 
four decades Nowadays, it is known that trans- 

port properties of conductors with extremely small cross 
sections, arc substantially controlled by the scattering 
of electrons at random inhomogencitics of the conduc- 
tor boundaries @-[n|. Many of transport properties of 
mesoscopic guiding systems can be modeled as surface- 
disordered waveguides. 

Recent numerical studies of quasi-one-dimensional 
(quasi-lD) surface-disordered systems fl2l - [l4j have re- 
vealed a key difference from the standard models with 
bulk random potentials. Specifically, it was found that 
transport properties of waveguides with surface scatter- 
ing essentially depend on multiple characteristic lengths, 
in contrast with the conductors with bulk scattering for 
which the single-parameter scaling [l5[ occurs. Accord- 
ing to the single-parameter scaling, all scattering prop- 
erties of finite conductors with bulk random potentials 
are determined by a single parameter: the ratio of the 
conductor size to the localization length found for the 
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corresponding systems of infinite size. However, the scat- 
tering properties of surface-disordered waveguides are de- 
termined by a non-isotropic character of scattering in the 
"channel space." 

For surface-corrugated waveguides with a large number 
of conducting channels Nd 3> 1, it was emphasized fl6l - 
l2lT | that the mode-attenuation length L n of the nth prop- 
agating mode (or mean free path of the nth conduct- 
ing channel) displays a rather strong dependence on the 
mode index n. Specifically, the larger the number n the 
smaller the corresponding mode-attenuation length and, 
as a consequence, the stronger is the scattering of this 
mode into the others. As a result, there emerges the 
hierarchy of mode-attenuation lengths [l6| : 



LN d < Ljsr d -\ < . . . < L 2 < L\. 



(1.1) 



The smallest length L^ d belongs to the highest (last) 
channel with the mode index n = N4 3> 1, while the 
largest length, L\, corresponds to the lowest (first) chan- 
nel with n = 1. Because of this hierarchy of lengths, 
a very important phenomenon can arise, that is the 
coexistence of ballistic, diffusive, and localized trans- 
port mm. For example, for a conductor of length 
L which is inside the set of lengths L n , 



L Nd < ... < L < 



(1.2) 



the lowest modes can be in ballistic regime, however, the 
intermediate and highest modes can exhibit the diffusive 
and localized behavior, respectively. 
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Accordingly, for quasi-lD surface-disordered waveg- 
uides one can state the following: (i) The average mode- 
transmittances (Tm), can be very different for different 
conducting channels, (ii) All propagating modes are 
mixed because of inter-mode transitions, therefore, the 
transmittance (Tm) of any given nth mode depends on 
the scattering into all modes, (iii) The total average 
transmittance, (T), of the waveguide contains the im- 
print of all the average mode-transmittances and, there- 
fore, can exhibit an coexistence transport in which the 
phenomenon of coexistence of different mode-transport 
regimes takes place. 

In this paper the above statements are quantitatively 
validated by means of two special models of two-mode 
waveguides with lateral symmetric rough boundaries 
(SRB) and antisymmetric rough boundaries (ASRB), see 
Fig. [1] The SRB waveguide represents a system in which 
the attenuation of the propagating modes is associated 
with the scattering within the same mode (intra-mode 
scattering). Comparatively, the attenuation in the ASRB 
waveguide arises due to the intra-mode scattering and 
due to the scattering of waves from one mode into the 
other (inter-mode scattering). 

Although main attention is paid to the conduction of 
electrons, the results are also applicable to the propaga- 
tion of classical waves. The latter problem has even a 
longer history, it naturally arises in the analysis of spec- 
tral and transport properties of optic fibers, acoustic and 
radio waveguides, remote sensing, shallow water waves, 
multilayered systems, photonic lattices and so forth [22I — 

n. 

This paper is structured as follows. In Section Hfl the 
geometry and general statistical properties of the surface- 
disordered waveguides are described. Here the Green's 
function of the waveguides and standard linear response 
formula are outlined, which are used to obtain the trans- 
mittance (or dimensionless conductance). Also, the for- 
mulas for mode-attenuation lengths of waveguides are 
discussed in detail. In Sections Mil and IIVI the trans- 
port properties of the SRB and ASRB waveguides with 
two propagating modes are analyzed. Here the expres- 
sions for the mode and total average transmittances are 
derived. The transport regimes of the total average 
transmittance are then analyzed; one being the coexis- 
tence regime, in which the coexistence of different mode- 
transports emerges. Finally, the conclusions are pre- 
sented in Section |Vl 



II. GENERAL STATEMENT OF THE PROBLEM 

A. Surface-disordered waveguides 

This study considers open surface-disordered waveg- 
uides (or conducting wires) of length L and average width 
d, stretched along the x- and z-axes, respectively; for 
quasi-lD geometry, naturally, L d. Two different con- 
figurations of the opposite boundaries result in the sym- 



metric and antisymmetric rough boundaries waveguides, 
see Fig. Q] The SRB waveguide has inhomogeneities that 
are symmetric with respect to its straight central guid- 
ing axis (dot-dashed line in Fig. [IJ. Thus, the inhomo- 
geneities give rise to the varying width of the SRB waveg- 
uide. In contrast, the width of the ASRB waveguide stays 
constant along the waveguide despite the rough bound- 
aries. Here the inhomogeneities are physically equivalent 
to a weak bending of waveguide. 

The lower and upper surfaces of the waveguides are de- 
scribed, respectively, by the rough boundaries z = o~£(x), 
and z = d^fcr^(x). The upper surface of the SRB (ASRB) 
waveguide is indicated by the minus (plus) sign of the 
symbol "=p" . Here a is the root-mean-square roughness 
height, which is assumed to be identical for both bound- 
aries. Hence, the waveguides occupy the region, 

-L/2< x < L/2, (2.1a) 
cr^(x) < z <d=F v£(x), (2-lb) 

of the (x, z)-plane. The random function £(a;) describes 
the roughness of the boundaries and is assumed to be sta- 
tistically homogeneous and isotropic, with the statistical 
properties of zero mean and unit variance, 

(£(*)) = 0, (2.2a) 
<£ 3 (*)> = 1, (2.2b) 
(Z(x)Z(x')) = TW(\x-x'\). (2.2c) 

Here the angular brackets represent the statistical aver- 
aging over different realizations of the surface profile £ (x) . 
The binary correlator W(x) decreases on the scale R with 
the normalization W(0) = 1, where R is the roughness 
correlation length. It should be reiterated that the mi- 
nus (plus) sign of the correlator corresponds to the SRB 
(ASRB) waveguide. Since lower and upper boundaries 
have the same rough profile, the binary correlator (|2.2cp 
describes the correlations within each boundary, as well 
as the cross-correlations between them. 



B. Transmittance and Green's function 

The transport properties of the disordered waveguide 
will be characterized by its transmittance or, equally, by 
its dimensionless (in units of e 2 /irh) conductance T(L). 
Within the standard linear response theory, as indicated 
by R. Kubo [27|, this quantity is expressed as 

4 r L / 2 r d 
T(L) = — p2 / dx dx' / dz dz' 

L J -L/2 Jo (2 3) 

dg(x,x';z,z') dQ*{x,x';z,z') 
dx dx 1 

where the asterisk "*" denotes complex conjugation. The 
retarded Green's function Q(x, x'; z, z') obeys the follow- 
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ing boundary- value problem [201 ] : 
d 2 d 2 \ 

— U(x, z) G{x, x'; z, z') = 5{x — x')S(z — z'), 

(2.4a) 

g(x,x';z = 0,z') = 0, (2.4b) 
£(x,x';z = eZ,z') = 0. (2.4c) 

Here <5(x) and S(z) are the Dirac delta- functions. The 
wave number k is equal to the Fermi wave number for the 
electrons within the isotropic Fermi-liquid model. For 
electromagnetic waves with both frequency uj and TE 
polarization, propagating through a waveguide with per- 
fectly conducting walls, k = lo/c. 

The equations (|2.4[) specify a Dirichlet boundary-value 
problem, however, with two remarkable features: the 
"bulk" scattering potential U{x,z), and the flat bound- 
aries of the waveguide at z — and z — d. These features 
emerge after the coordinate transformation of the initial 
Dirichlet boundary-value problem, in which there is no 
"bulk" potential while the scattering is only caused by 
the roughness of the boundaries. The idea of this method 
was first discussed by Migdal (28|, and has since been 
frequently used in theories of wave or electron surface- 
scattering. Details of the application of this method 
for "flattening" the rough boundaries, can be found in 
Ref. HtJ The form of the potentials U(x, z) for the spe- 
cific waveguides used in this study can be extracted from 
the aforementioned reference. It should be stressed that 
this potential can not be treated as a completely random 
potential and modeled, for example, by random matrices. 
A close inspection of this potential shows that random 
function £(x) enters the potential in a rather complicated 
way, however, by no means, strong correlations between 
matrix elements can not be neglected. 

For a further discussion we have to indicate that an ap- 
proximate expression of the potential for the SRB waveg- 



uide contains several terms that can be joined together in 
three groups. These groups depend on cr£ (x) /d, at;' (x) / d, 
and <j 2 £' 2 {x)/d 2 , 



U(x,z) w U 



d ' d ' d 2 



-, z 



(2.5) 



An exact expression for the ASRB waveguide contains 
other terms, which, however, can be arranged in the 
groups that depend on a£'(x)/d and a 2 £' 2 (x) / d 2 , 



U(x,z) = U 



aC(x) a 2 C 2 {x) 
d ' d 2 ' 



(2.6) 



Here the prime to the function £(x) denotes a derivative 
with respect to x. 

As is typical, after the application of the Green's theo- 
rem, the problem can be reformulated in the form of the 
Dyson equation. In order to perform the averaging of the 
Green's function, one can apply one of the well-known 
perturbative methods. For instance, the diagrammatic 
approach developed for surface-disordered systems [22j], 
or the technique developed in Ref. may be equally 
utilized. Both methods take adequately into account the 
multiple scattering from rough boundaries and allow one 
to formulate a consistent perturbative approach with re- 
spect to the scattering potential. After straightforward 
but cumbersome calculations (see Ref. [TH for details) , the 
following average Green's function is obtained, 



(g(\x - x'\;z,z')) =^sin (—j-) sin 



exp[(ifc„ - l/2L n )\x - x'\] 



(2.7) 



Here Nd = {kd/ir^ is the total number of the propagating 
normal modes in the flat waveguide, which is determined 
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by the integer part [. ..] of the mode parameter kd/n. 
Due to the quantization of the transverse wave number 
k zn = im/d, the quantum value of the longitudinal wave 
number is given by 



k x „ = k„ = ^/k 2 - (nn/d) 2 , n = 1, 2, 2V d . 

(2.8) 

All other waveguide modes with n > Nd are evanescent 
with imaginary values of k n , and do not contribute to the 
transport properties. The average Green's function (|2.7[) 
differs from the unperturbed one, in the appearance of 
the mode- attenuation length L n . 



C. Mode-attenuation length 

The mode-attenuation length (or electron total mean 
free path) L n describes the scattering from the nth prop- 
agating mode into the others. From general theory of 
quasi-lD scattering systems it follows that the mode- 
attenuation length is determined by the backward Ln , 
and forward L„ scattering lengths [22j , 

T~ = T(b) + TU)' ( 2 ' 9 '' 

L n L n ' L n ' 

In accordance with Ref. [2(| both backward and forward 
scattering lengths consist of two fundamentally different 
partial lengths. One of these partial lengths is related to 
the amplitude and gradient mechanisms of surface scat- 
tering, whereas the other partial length is associated with 
the square- gradient mechanism (the connection between 
these partial lengths and the scattering potential U(x,z) 
is discussed below). Therefore, backward and forward 
scattering lengths are given by 



1 

i-in 



■(f) 



1 



r(A,G)(b) 
-Lin 



r (SG)(6)' 
-Lin 



1 



1 



(A,G)(f) 



-(«?)(/)' 



(2.10a) 
(2.10b) 



The partial lengths associated with the backscattering 
read 



1 



(A,G)(6) 



L 



a- 
7> 



1 



L 



(SG)(b) 



N d 

E 

n'=l 
4 N d 

-T 

d 4 ^ 

n' = l 



A, 



-W(k n + k n r), (2.11a) 



Br, 



kn k r 



-S(k n + k n >). (2.11b) 



The partial lengths associated with the forward scatter- 
ing are similar to the backward ones but with the func- 
tions W and S depending on the argument k n — k n ', i.e., 



a 2 A , 

V v ^ /Inn/ 



£ YY^W{k n -k n/ ), (2.12a) 



L (A >G )(/) ' d 6 knkn 

1 er 4 Nd B 

AsgW) = kjt> S{ - kn ~ kn,) ' (2 ' m) 

-Lin n ' = l 



Factors A nn > and B nn ' depend upon the type of sym- 
metry exhibited by the opposite rough boundaries of the 
waveguide. Below, they are discussed and, in Table U 
their values are explicitly shown for the SRB and ASRB 
two-mode waveguide. 

The equations (|2.11[) and (|2.12j) contain two differ- 
ent spectral functions. First, the roughness-height power 
spectrum is included in Eqs. (|2.11a[) and (|2.12a[) . This 
function is defined as the Fourier transform of the 
roughness-height binary correlator W(x), see Eq. (|2.2cj) ). 



W{k x 



dx exp(— ik x x) W(x). 



(2.13) 



Second, the formulas in Eqs. (|2.11bp and ()2.12b[) contain 
the so-called roughness-square-gradient power spectrum, 



S(k x ) 



dx exp (— ik x x) W" 2 (x) 



(2.14) 



where the double prime to the function W(x) denotes 
a second derivative with respect to x. Since W{x) and 
W" 2 (x) are real and even functions of x, their Fourier 
transforms are even and real functions of the longitu- 
dinal wave number k x . It should be also stressed that 
according to rigorous mathematical theorem, the power 
spectra are a non-negative functions of k x for any real 
random process £(x). 

The spectra (|2.13|) and (|2.14|) arise when deriving the 
correlator of the scattering potential U{x,z) in the k x 
representation. The correlator of U (x, z) emerges when 
averaging the Green's function within the perturvative 
approach [l8|]. Specifically, the terms in U(x,z) that de- 
pend upon the roughness amplitude £(x) and the rough- 
ness gradient £'(x) give rise to the terms in the correla- 
tor associated with W(k x ), see the structure of U(x,z) 
in Eqs. (|2.5[) and (|2 . 6[) . The terms that depend upon 
the roughness square gradient £,' 2 (x) lead to the terms 
in the correlator associated with the Fourier transform 
of (C(x)C(x')} 2 = W 2 (x), which in fact is S(k x ). It 
should be stressed that through the integration by parts 
the power spectrum of the roughness gradient £'(x) can 
be reduced to W(k x ). However, it is not possible to do 
the same for the power spectrum S(k x ). This fact reflects 
a nontrivial role of the terms in the scattering potential 
containing £,' 2 (x). Therefore these terms are associated 
with a specific and independent square-gradient surface- 
scattering mechanism. 

Accordingly, the inverse scattering lengths given in 
()2.10j) arc expressed as a sum of two terms that describe 
scattering governed by different surface scattering mech- 
anisms. Thus, the amp/i£uc?e-scattcring and gradient- 
scattering mechanisms are associated with the terms 

1 / L (A, G)(b) ^ 1 / L (A,G)(f) giyen in Eqf ^ (J^Y^ and 

(|2.12ap . while the sgMare-gradieni-scattering mechanism 

is associated with the terms l/Ll SG '"- b ' > and l/ii SG ' ) ^' > 
given by Eqs. (|2.11b[) and (|2.12b|) i. 

It should be also stressed that, as functions of the cor- 



relation length R, the mode-attenuation lengths L 



(A,G)(b) 
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and L^' G ^^' can behave very different from L^f^^ and 
£^(SG)(/) j-, ecause f their dependence upon the different 
power spectra W(k x ) and S(k x ). Owing to those peculiar 
behaviors, the competition between surface-scattering 
mechanisms emerges and, for some range of values of 
R, the product a 4 S(k x ) can greatly increase its value in 
comparison with the quantity a 2 W{k x ), even for small a. 

Thus, the lengths L^ SG ^ b ^ and i}^ can be compared, 



(A,G)(b) 



and 



and even can be smaller than the lengths L 
^CA,G)(/) detailed discussion about this competition 
between surface-scattering mechanisms is presented in 
Section inTl for the SRB waveguide, and in Section[TVl for 
ASRB waveguide, with the specific power spectra (|2.22l) , 
which are associated with the Gaussian correlator. 

The approximation by which one arrives to Eqs. (|2.7p . 
(|2.9p - (|2.12[) . and the resulting domain of applicability 
of those equations were discussed in Ref. [TH Specifically, 
it was shown that the validity of the results is restricted 
by two independent criteria of weak surface scattering, 



A„ < 2L n , 
R < 2L n . 



(2.15a) 
(2.15b) 



Here the cycle length, A n , is the distance between two 
successive reflections of the nth mode from the rough 
boundaries, 



A„ = k n d/(irn/d). 



(2.16) 



The criteria expressed in Eq. (|2.15[) imply that the waves 
are weakly attenuated over both the correlation length R 
and cycle length A„. Additionally, they restrict the cor- 
rugations to be small in height, a <C d. It should be also 
noted that statistical averaging is reasonable if the cor- 
relation length is much less than the waveguide length, 
fi<L. Therefore, the length L must obey requirements 
similar to Eqs. (|2.15[) formulated with respect to L n . All 
these limitations are common in surface scattering the- 
ories that are based on an appropriate perturbative ap- 
proach, see for example, Ref. 1221 



D. Reduction to two-mode waveguides 

In what follows, this study is restricted to waveguides 
with two propagating modes, 



N d = lkd/n} = 2. 



(2.17) 



In this case the mode parameter is confined within the 
interval, 

2 < kd/n < 3, (2.18) 
and the longitudinal wave numbers are given by 



ki = yjk 2 - (ir/d) 2 , k 2 = ^Jk 2 - (2n/d) 2 . (2.19) 



Then, the factors A nn i and B nn i, in Eqs. (|2.11[) and 
()2.12j) . which are drawn from Ref. [2(| are displayed in 



TABLE I. Matrices of factors A nn i and B nn i for the symmet- 
ric and antisymmetric two-mode surface-disordered waveg- 
uides. 



SRB 



ASRB 



An An 

A%\ A22 

Bu B12 

B21 B22 



4tt 4 
64tt 4 



16tt 4 



16tt 4 







18 






(3+4tt 2 ) 2 
18 



71-72 

8tt 4 



Table HI The diagonal and off-diagonal elements are as- 
sociated, respectively, with intra-mode scattering (scat- 
tering within the same mode) and inter-mode scattering 
(scattering from one mode into the other). Also, one can 
realize that the diagonal and off-diagonal elements of the 
matrices A nn i are related to the amplitude-scattering and 
gradient-scattering mechanisms, respectively. 

Bearing in mind this physical meaning of the entries 
A nn i and B nn i , the structure of matrices in Table U re- 
veals remarkable differences between the SRB and ASRB 
waveguides: 

(i) The attenuation in the SRB waveguide is deter- 
mined by the competition between the amplitude- 
scattering and square-gradicnt-scattering mecha- 
nisms, which are both associated with the intra- 
mode scattering (in two-mode waveguide). There- 
fore, the propagating modes in the SRB model 
should be regarded as independent, and, conse- 
quently, the SRB two-mode waveguide represents 
two independent lD-surfacc-disordcred channels. 
As is known, the transport properties of ID dis- 
ordered systems are described by the theory of An- 
derson localization (see, e.g., Ref. Il6l and references 
therein). This theory has been already employed 
for the analysis of SRB single- mode waveguide [3(| ■ 
Now, the theory of ID Anderson localization should 
be adapted for describing the transmittance of two- 
mode waveguide. According to this theory, the 
transport in each channel is completely specified 
by the scaling parameter that is the ratio between 
the sample length and corresponding backscatter- 
ing length. Remarkably, the transmittance of any 
ID disordered structure does not depend on the 
forward scattering. This approach is presented in 
Section [TTTl where it is showed that L/Ln^ becomes 
the scaling parameter of the transport through the 
nth channel. 

(ii) Otherwise, the attenuation in the ASRB waveg- 
uide is determined by the interplay between the 
gradient-scattering and square-gradicnt-scattering 
mechanisms. Since the former is related to the 
inter-mode scattering, the channels can become 
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mixed and the theory that accounts for the trans- 
port properties of the ASRB waveguide should dif- 
fer from that used in SRB waveguide. A general 
theory for the transmittance of waveguides that 
considers the interplay between both mechanisms 
is not presented here. However, the theoretical ap- 
proach developed in Refs. [TH and [2^ allows one to 
predict many of important characteristics of trans- 
port in different regimes. Here this approach is 
applied to the ASRB waveguide with two chan- 
nels. As a result, it becomes possible to predict the 
interval of parameters in which only one of these 
mechanisms prevails. On the one hand, when the 
intra-mode scattering governs the attenuation, one 
deals with the same situation that was found in 
the SRB waveguide and, therefore, the transport 
properties may be specified by the scaling parame- 
ter that is the ratio between the sample length and 
corresponding backscattering length. On the other 
hand, when the inter-mode scattering determines 
the attenuation the transmittance of the nth prop- 
agating mode is specified by the scaling parameter 
L/L n , which involves both backward and for- 
ward L)i scattering lengths. The approach that 
considers intra-mode and inter-mode scattering in 
the ASRB waveguide is developed in Section HVl 

The explicit expressions for the mode backward and 
forward scattering lengths, of our two-mode waveguides, 
are useful for continuous reference. From Eqs. (|2.10p - 
(|2.11[) , with the factors A nn i and B nn i from the column 
SRB in Table U the necessary backscattering lengths are 
given by 

' 4.^2^+9.2 (2.20a) 
1$) d 6 k\ d 4 k\ y ' 



1 

lb) 



d e 



k 2 



d A k\ 



(2.20b) 

Also from Eqs. (|2TTU|) - (j2~T2"j) . but with the factors A nn > 
and B nn i corresponding to the ASRB waveguide, the ex- 
plicit expressions of the lengths read, 



L 



1 



if) 



= 167T 



= 16tt 



a 2 W{k 1 + k 2 ) 7T 4 a 4 S(2fci) 
d 6 kik 2 + Yd 4 k\ ' 



(2.21a) 



a 2 W{k 1 -k 2 ) 7T 4 a 4 5(0) 



(2.21b) 



L 



CO 



1fi 4 a 2 W(k! + fc 2 ) , p i cr i S{2k 2 ) 

107T —^r ; — ; V OTT 



d e kik 2 



d 4 k\ ' 



(2.21c) 



L 



(f) 



d 6 kik 2 



d 4 k 2 



(2.21d) 



In this paper we consider the situation for which the 
correlator W(x) can be approximated by the Gaussian 
form, W(x) = cxp (— x 2 /2R 2 ). After the substitution 
of this correlator into formulas (|2.13|) and (|2.14p . for 
any degree of roughness, the power spectra appearing 
in Eqs. ([2301) and P^T|> arc given by 



(2.22a) 



W(k x ) = V2irRexp[-(k x Ry /2], 



S(k x ) 



16R 3 



[{k x Rf - 4(k x R) 2 + 12] 



x exp[-(fc x i?) 2 /4]. 



(2.22b) 



These spectra have their maximum at k x 
W(0) - R and S(0) - R~ 3 . 



0, but with 



III. TWO-MODE SRB WAVEGUIDE 



A. SRB: Competition between surface-scattering 
mechanisms 



In order to discuss the competition between various 
mechanisms of surface scattering, we should distinguish 
between two regions for the correlation length R. These 
regions are denoted as the region of small-scale rough- 
ness {kR 1) and the region of large-scale roughness 
(kR ^> 1). Both regions correspond to weak correlations 
between successive reflections of the wave from rough 
boundaries (R <C A n ). Note that due to the last re- 
quirement the second of the weak-scattering conditions 
(|2.15p is satisfied automatically when the first one is met, 
i?< A„ < 2L n ,L. 

Small-scale-roughness - In this region the surface 
roughness can be regarded as a delta-correlated random 
process of the white-noise type. Taking into account the 
evident relationship 1 < fcA„ and the weak-scattering 
conditions (|2.15[) . one can get the following inequalities 
to specify this region together with requirements of ap- 
plicability of the theory: 



kR < 1 < kA n < 2kL n . 



(3.1) 



Under the condition (|3.1[) . the argument of the power 
spectra (|2.22[) is much less than the scale of their de- 
crease, k x <C R ■ Therefore, one can write 



(3.2a) 
(3.2b) 



W(k x ) w W(0) = V2ttR, 
S(k x ) « 5(0) = 3^/4-R 3 . 



In this case, when R decreases, the spectrum S(k x ) 
increases as fast as R~ 3 , whereas W(k x ) decreases as 
R. As one can see from Eqs. (|2.20[) . in spite of the fact 
that the first term is proportional to a 2 while the sec- 
ond is proportional to cr 4 , for any value of the roughness 
height cr, there is a region of small values of R, where 
the square-gradient-scattering mechanism predominates 
(provided the condition (|3.1[) is fulfilled). Specifically, if 
the correlation length R is smaller than some crossing 
point R cn , the first term in Eqs. (|2.20ap and ()2.20bP can 
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be neglected and the mode-backscattering lengths are ap- 
proximated to 



-4rv « 12.2 fci for fci? < fci? ci < 1, 



(6) 



(3.3a) 



(6) 



(a/d) 4 

133.3 fc 2 ;, for fci? < fci? c2 < 1. 

(fc 2 i?) J 



(3.3b) 

Otherwise, when i? is larger than i? cn , the square- 
gradient (second) term in Eqs. (|2.20a|) and (|2.20b|) can 
now be neglected and the mode-backscattering lengths 
are approximated to 



Li 



L 



(b) 



10.0 fci 



160.4 h 



(a/d) 2 (hR) 
(kid/w) 4 

for kR cl < kR < 1, (3.4a) 

(q/rf) 2 (fcgg) 

(fc 2 d/7r) 4 

for kR c2 < fci? < 1. (3.4b) 

The dimensionless crossing point fci? cn of the nth 
backscattering length can be located either on the 
border between the regions of small- and large-scale- 
roughness, or inside the first region. In the former case, 
kR cn ~ 1, and within whole region ()3.1[) . the mode- 
backscattering lengths are contributed mainly by square- 
gradient-scattering, i.e. by Eq. (|3.3p . In the latter case, 
fci? c i and fci? C 2 are found by searching for the intersec- 
tion of the asymptote (|3.3a|) with (|3.4a[) , and (j3.3b[) with 
(I3.4bp . Thus, the crossing points read 



fci? c i ~ fci? c 



0.3 kd 



(3.5) 



Large-scale-roughness - This region arises when 
the correlation length i? becomes much larger than the 
wave length 27r/fc, but still remains much less that the 
cycle length A„, 



1 < fci? < fcA„ < 2kL r 



(3.6) 



Here the square-gradient term in Eqs. (|2.20a[) and (12.20b)) 
can be neglected and the mode-backscattering lengths are 
approximated to 

(a/d) 2 (hR) 



1 

lb) 



10.0 fci 



160.4 fc 2 



(M/tt) 4 

for 1 < fci? < fcAi, 
(a/d) 2 (k 2 R) r 

77 J] , 4 CX P - 

(fc 2 d/7r) 4 
for 1 < fci? < fcA 2 . 



exp [-2(fcii?) 2 ] 

(3.7a) 
2(fc 2 i?) 2 ] 

(3.7b) 



B. SRB: Total and mode transmittances 

Based on the diagonal form of matrices A nn i and B nn > 
of the SRB waveguide, this study argues that the at- 



tenuation of waves arises due to the intra-mode scat- 
tering only (see Table. HJ). Consequently, the two-mode 
waveguide can be treated as two independent lD-surface- 
disordered wires with their corresponding backscattering 

lengths and L\ . Thus, the known theory of ID lo- 
calization becomes pertinent here. In adapting this the- 
ory to the present situation, it should be realized that the 
scaling parameters L/L^ and L/L 2 b ^ determine, respec- 
tively, the average transmittance of the first and second 
mode. Specifically, the transmittance of the nth mode 
(n = 1,2) can be taken, for example, from Ref. [l6|. It 
reads 



(Tm(L/L^)) 



1 / L 



: 2y^ V 



41 



(b) 



-3/2 



exp 



L 



41 



('■') 



cosh z 



exp 



L 



The total transmittance is given by, 



(T(L))=Y J {T M (L/L^)). 

n=l 



(3.8) 



(3.9) 



Eqs. (|3.8j) and (|3.9[) can be directly found from the 
general expression (|2.3[) with the use of well developed 
methods, such as, e.g., the perturbative diagrammatic 
technique of Berezinski [2, EJ , the invariant imbedding 
method [HI, [33| or the two-scale approach [3(1 H3, HH . 

Here it should be stressed that the ratio L/Ln^ in- 
cludes the contribution of both the amplitude-scattering 
mechanism and the square-gradicnt-scattering mecha- 
nism through the first and second terms in Eqs. ()2.20a[) 
and (|2.20b|) . 

The mode transmittance (|3.8j) exhibits the ballistic be- 
havior for large backscattering length, 

(T m (L/lW))*1-L/lW for L«4 6 >. (3.10) 

In this case the nth conducting channel is practically 
transparent. On the contrary, the the mode transmit- 
tance displays exponential decrease as the waveguide 
length L exceeds 4i£ 



-(b) 

J n , 



(Tm(L/L^)) 



7T 5 / 2 / L 



16 



-3/2 



for 



exp 



L 



L^«L. (3.11) 



The dependence of the function (|3.8j) and its asymptotes 

(|3.10[) and (|3.1ip on the scaling parameter L/L^n are 
shown in Fig. [2] 

In accordance with two mode-transport regimes 
p.lOp . (|3.11l) and due to the hierarchy of the mode- 
backscattering lengths, L\' < (this hierarchy can 
be seen when comparing Eq. (|2.20a[) with (|2.20bl) . or 
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when the smallest mode-backscattering length L 2 
is larger than L. As a consequence, both conduct- 
ing channels are open having almost unit transmit- 
tances (|3.10p . The waveguide is almost perfectly 
transparent. Its total transmittance (|3.9[) approxi- 
mately equals to the total number of the propagat- 
ing modes, 

(T(L)) » 2. (3.17) 



FIG. 2. (Color online) SRB waveguide: the mode transmit- 
tance (|3.8[1 (solid line) and its asymptotes (|3.10p and (|3.11[) 
for ballistic and localized regimes, respectively (dashed lines). 



more directly through the couples of Eqs. (|3.3[) . (|3.4p 
and (|3.7[) ). the total transmittance (|3.9j) has the follow- 
ing three transport regimes: 

(i) In the regime of localization when the largest 
mode-backscattering length is smaller than the 
waveguide length L, 



1<L/L[ b) <L/Li h \ 



(3.12) 



both propagating modes are strongly localized, and 
their transmittances are exponentially small. The 
total transmittance (|3.9[) is approximately equal to 
the transmittance of the first mode obeying the 
asymptote p. lip . 



(T(L)) « (T M (L/L^)) 



r 5/2 



L 



-3/2 



16 \ALf ] J 



4L 



(3.13) 

The waveguide is non-transparent in this regime. 

(ii) The coexistence regime arises when the smallest 
backscattering length L 2 is smaller, while the 
largest backscattering length is larger than the 
waveguide length L, 



L/L[ b) < 1 < L/L {b) . 



(3.14) 



In this case, the first mode manifests the ballistic 
behavior (|3.10p . while the second mode is still lo- 
calized in line with Eq. p. lip . Therefore, as before, 
the total transmittance (|3.9p is determined by the 
transmittance of the first mode, however governed 
by the ballistic asymptote (|3.10p . 

(T(L)) « (T M (L/L®)) « 1. (3.15) 

(iii) The ballistic regime emerges under the conditions 

L/L[ b) < Ljhf < 1, (3.16) 



C. SRB: Transmittance vs correlation length 

Here we discuss the dependence of the mode- 
backscattering lengths and transmittances on the dimen- 
sionless correlation parameter kR. For this, we plot the 
ratio L^' /L 2 > as we ^ as ^ e scann S parameters L/L] 1 ' 
and L/L 2 \ see Figs. [3] (a) and (b). These plots are com- 
puted from the expressions for iJh 1 given in Eqs. (|2.20a|) 
and (|2.20bp with power spectra (|2.22p . If R decreases, 
the lengths L n decreases as R 3 whereas A ra remains fixed, 
thus, the smallest value of R is restricted by the criterion 
(|2.15a|) : this criterion is fulfilled since A\/2Li w 0.01 <C 1 
and A2/2L2 ~ 0.04 <§C 1 at the left boundary of the plots, 
kR = 0.01. Because the plots illustrate the case of weak 
correlations (i{ < A n ), the condition (|2. 15h>[) is satisfied 
automatically if (|2.15ap is met. 

Figure [3] (a) exhibits the crossover from the square- 
gradient to amplitude scattering. Specifically, with an 
increase of kR the ratio / starts with the value, 



L 



(6) 
1 

(b) 



10.9 



for kR < kR c2 , (3.18) 



and grows in accordance with the relation, 



(b) 
l 

(6) 



16 



fei 



exp [6tt 2 (R/c1) 2 ] 

for kR cl < kR < fcA 2 . (3.19) 



The asymptotic behaviors (|3 . 18[) and (|3 . 19[) are obtained, 
respectively, from Eqs. (|3.3p and (|3.7p . Note that the 



hierarchy > L 2 remains for any value of the pa- 
rameters. With the curves in Fig. [3] (b), the different 
transport regimes defined in Eqs. (|3.12p - (|3.17p can be 
identified. 

The mode transmittances (|3.8p and total transmit- 
tance p.9p are plotted in Fig. [3] (c). Because of the 
parameters used in the plots, the boundary between the 
localization and the coexistence regime is at kR ps 0.02. 
With further increment of kR the propagation regime 
becomes ballistic after the point kR w 0.06. Figure [3J 
(c) clearly demonstrates that, due to the hierarchy of 
backscattering lengths, L 2 < L>^\ the transmittance of 
the second mode is always smaller than the transmittance 
of the first one, 



(T M (L/L 2 b) )) < (T M (L/L[ b) )). 



(3.20) 
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3. (Color online) SRB waveguide: Mode-backscattering 

and transmittances, (TM n ) — (Tm(L/ L 1 ^)), 
(T(L)), vs dimensionless correlation length, kR, for o/d = 
0.0025, fcd/vr = 2.54, (ka « 0.02), and L/d = 1600. The 
crossing points fcJ? c i ~ 0.13 and fcJ? C 2 ~ 0.12 are indicated 
on the top of the upper plot. The boundary between the 
localization and the coexistence regime is at kR « 0.02 and 
with further increment of kR the propagation regime becomes 
ballistic after the point kR « 0.06. 



Therefore, within the localization (|3. 12|) and coexistence 
(|3.14p regions the total transmittancc (|3.9I) is mainly con- 
tributed by the first term, i.e. by the first-mode trans- 
mittance, 



(T(L)) « (T M (L/L^)) 



(3.21) 



Within the localization regime, (T(L)) behaves as in- 
dicated in Eq. (|3.13[) . Within the coexistence regime 
the total transmittance is described by Eq. (|3.15p . For 
the ballistic regime p,16p . the approximation given by 
Eq. p.lOp describes well the behavior of the mode trans- 
mittances. Correspondingly, the value of the total trans- 
mittance tends to the total number = 2 of the prop- 
agating modes, see Eq. (|3.17[) . 

Thus, with the variation of roughness correlations one 
can realize all three regimes, which is inherent for the 
transport through the two-mode SRB waveguide. Note 
one more that coexistence regime (|3.14[) , (|3.15[) arises 
only due to the fact that in the SRB waveguide the hierar- 
chy of the mode-backscattering lengths is always present. 



D. SRB: Transmittance vs mode parameter 

From the experimental viewpoint, it could be more fea- 
sible to study the transport properties of a waveguide by 
fixing its geometrical parameters such as length, width, 
etc., whereas the wave number k is varied within some in- 
terval. In Fig. [4] we present three pairs of frames formed 
by plots of mode-backscattering lengths and transmit- 
tances as functions of the mode parameter kd/n (di- 
mensionless wave number) within the range of first two 
modes. The pairs (a)-(b), (c)-(d) and (e)-(f) mainly illus- 
trate, respectively, the localized p,12p , coexistence (|3.14p 
and ballistic p,16p transport regimes. Here the different 
regimes are illustrated by selecting three specific values of 
the correlation length R whereas the length of the waveg- 
uide L stays constant; when kd/n = 2.54, the mode- 
backscattering lengths and transmittances corresponds 
to those presented in Fig. [3] at the points kR = 0.014, 
kR = 0.04 and kR = 0.3. In this way, inside each regime, 
the relative influence of the scattering mechanisms is dif- 
ferent. Another possibility to arrive to different transport 
regimes, which is not shown here, would consist in select- 
ing only one value for kR and three different values of the 
length L. 

Figure U] (a), (c) and (e) shows the plots of Ln^/L 
computed from Eqs. (|2.20a[) and (|2.20b[) . as well as 



L 



A,G)(b) 



/L, which also comes from these equations but 
when the square-gradient term is neglected. If kd/ir in- 
creases, the nth attenuation length also increases from its 
value equal to zero at kd/ir = n, where the nth normal 
mode opens (n = 1,2). Note that our results, which are 
obtained within the second order approximation in the 
perturbation potential, can only indicate the position of 
these zeros and cannot correctly describe their line shape 
in the narrow vicinity of integer values of kd/ir. The 
reason is that at these resonance points the longitudinal 
wave number k n vanishes and, as a consequence, the cor- 
responding inverse attenuation length diverges. Conven- 
tionally, this problem is solved within the self-consistent 
approach providing finite resonance values for 1/Ln (see, 
e.g., Ref. (22|). It should be mentioned that the curves 
associated with the first propagating mode do not man- 
ifest any resonance at kd/ir — 2, even with the open- 
ing of the second propagating mode. This occurs since 
both mechanisms that govern the scattering in two-mode 
SRB waveguide are associated with intra-mode transi- 
tions only. Remarkably, in frames (a) and (c) the dimen- 
sionless lengths L^'^^'/L and /L are much 
larger than the total lengths, L^/L and l\ b '/L. This 
is because the main contribution is due to the square- 



gradient-scattering lengths, L\ ^ >(u ' /L and L^"^^ ' / L. 



\ SG){b) /L and 4 S G)(b) / 
Therefore, for values of the mode parameter far from 
the resonances, the approximate behavior of Ln^/L is 
given by Eqs. p.3p (not plotted in Figure). Otherwise, 
in frame (e) , in which the main contribution is due to the 



amplitudc-gradicnt-scattering lengths, L 



(A,G)(b) 



/L and 
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L^'^ /L, the approximate behavior of Ln^ /L is given 
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FIG. 4. (Color online) SRB waveguide: Mode-backscattering 
lengths and transmittances versus mode parameter for a/d = 
0.0025, L/d = 1600, with R/d w 0.0018 for (a)-(b), w 
0.005 for (c)-(d) and R/d « 0.038 for (e)-(f). Note that 
for fixed values of the roughness height a and the mode pa- 
rameter within the interval 2 < kd/n < 3, fccr ranges the 
interval 0.016 ka 0.024; similarly, for fixed values of 
the correlation length R and the mode parameter within 
the interval 2 < kd/n < 3, the dimensionless correlation 
length kR ranges the interval 0.011 < kR < 0.017 in (a)- 
(b), 0.031 < kR< 0.047 in (c)-(d) and 0.23 < kR < 0.35 in 
(e)-(f). 



by Eqs. J53J). 

The mode transmittances (|3.8[) . and total transmit- 
tance (|3"1?|) arc plotted in Fig. H (b), (d) and (f). In 
frame (b) the first and second propagating modes are in 
localized transport regime; this fact is indicated by the 
small mode transmittances given by Eq. p. lip . The to- 
tal transmittance is approximately equal to the transmit- 
tance of the first mode, thus, it is given by Eq. f|3.13|) . As 
one can see in Fig. [3j with the increment of kR the mode- 
attenuation lengths also increases resulting in the condi- 
tion (|3.14p that is illustrated in frame (d). Here, the first 
and second propagating modes are, respectively, in ballis- 
tic and localized transport regimes. The coexistence situ- 
ation for the total transmittance is graphically indicated 
by a very small deviation of (T(L)) from (Tm(L/L{')}, 
when kd/n approaches 3. This behavior is marked by 
the ballistic asymptote (|3~T0|) for (T M (L/ Xf')) and by 

increasing rate of the curve (Tm (L / L^' ')) with the lo- 
calization asymptote (|3.11[) . The predominating contri- 
bution of the square-gradient-scattering mechanism can 
be noted visually in frame (b) and (d) by comparing the 
curve of (T(L)) with the curve (T^'°\L)), 



(T^ G )(L)) = ^(T M (i/L 



(A,G)(b) 



))• 



(3.22) 



71 = 1 



that ignores this contribution. In Fig. 0] (f) the first 
and second propagating modes arc in ballistic transport 
regime. When increasing kd/ir, this is portrayed by the 
sharp increase of the transmittance and, then, by the bal- 
listic behavior given by the asymptote (|3.10p . This step- 
like behavior indicates the openness of the channels with 
almost unit transmittance and, consequently, the total 
transmittance given by Eq. (|3.17p . as the total number 
of propagating modes. Since in frame (f) the main con- 
tribution is due to the amplitude-scattering mechanism, 
the behavior of the total transmittance follows the ap- 
proximation (|3.22p (see the overlapping of curves related 
to Eqs. (EH) and (jX2"2"]) ). 



IV. TWO-MODE ASRB WAVEGUIDE 

In contrast with the SRB waveguide, the theory for the 
ASRB waveguide can not describe the whole transition 
from ballistic to localized transport, which is expected 
for large enough L (numerical evidence of this transi- 
tion can be found, for instance, in Refs. [ID and HH). 
However, the analysis of the mode-attenuation lengths 
that is performed here for two mode waveguides, predicts 
many of important characteristics of transport. Thus, 
as is remarked in Sec. HI the attenuation in the ASRB 
waveguide is determined by the interplay between the 
gradient-scattering and square-gradient-scattcring mech- 
anisms, i.e., it is determined by the interplay between 
inter-mode and intra-mode scattering. Below, firstly, a 
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more detailed analysis of mode-attenuation lengths allow 
us to compute the interval of parameters in which only 
one of these mechanisms prevails, as well as their be- 
havior inside that interval. Secondly, the theory already 
developed for the SRB waveguide, in which the trans- 
port properties of the nth channel are specified by the 
scaling parameter L/Ln \ is applied when only the intra- 
mode scattering governs the attenuation and the channels 
can be considered as independent. Thirdly, the theory 
that accounts for the transport properties of the ASRB 
waveguide, when the inter-mode scattering determines 
the attenuation (when the channels become mixed), is 
developed. Finally, a small region of transition, where 
no analytical results exist, is only indicated. 



ASRB: Competition between surface-scattering 
mechanisms 



The competition between surface-scattering mecha- 
nisms is again discussed within the two intervals of cor- 
relation length _R, which are denoted as the region of 
small-scale roughness (kR <C 1) and the region of large- 
scale roughness (kR >• 1). Both regions correspond to 
weak correlations between successive reflections of the 
wave from rough boundaries (R -C A n ). See details about 
the definition of the regions of correlation length in sub- 
section MI Al 

Small-scale-roughness region - In this region, see 
Eq. (|3.ip . if the correlation length R is smaller than the 
crossing point R cn , the first term in Eqs. (|2.2ip can be 
neglected. Then, after the substitution of Eq. (|3.2b[) into 

the second term of Eqs. (|2.2ip one arrives to « 
l/Ln , and the following asymptotic formulas for the 
backscattering lengths are obtained 



lengths (|2.9[) are approximated to 



1 

Li 



— w 80.2 k 2 



(a/dfjhR) 
(kxd/n) 2 (k 2 d/Tr) 2 
for kR cl < kR < 1, 



1 



1 

Li 



— w — for kR c2 



(4.2a) 
(4.2b) 



The dimensionless crossing point kR cn of the nth 
length can be located either on the border between the 
regions of small- and large-scale-roughness, or inside the 
first region. In the former case, kR cn ~ 1, and within the 
whole small-scale roughness region the mode-attenuation 
lengths are described only by square-gradient-scattering 
contribution, i.e., by Eq. (|4.ip . In the latter case, kR c i 
and kR C 2 are computed by searching for the intersection 
of the asymptotes (|4.ip with (|4.2p . Thus, the crossing 
points read 



kR cl 



kR c 



0.36 



0.72 




(4.3a) 



(4.3b) 



Large-scale-roughness region - In this region, see 
Eq. (|3.6p . the square-gradient terms in Eqs. f|2.21[) can 
be neglected and the mode-attenuation lengths (|2.9p are 
approximated to 



1 

~L~x 



— w 40.1 k 2 



exp 



(a/d) 2 (kxR) 



(k 1 d/n) 2 (k 2 d/TT) 2 
; i 



(fci + k 2 ) 2 R 2 



+ exp 



(A* - k 2 ) 2 R 2 



for 1 < kR. < fcAi, 



(4.4a) 



w — , forl«M« kA 2 . (4.4b) 
L 2 L\ 



L 



(6) 



64.7 k x 



-4rr ~ 1035.9 fc 2 



{hxRf 
for kR < kR c x < 1, 
(a/dr 



(k 2 Rf 
for kR < kR c2 < 1. 



(4.1a) 



(4.1b) 



It is necessary to be aware that the present analysis of 
the competition is valid for the values of mode param- 
eter kd/w not very close to the resonances kd/ir = 1,2. 
Note particularly, that the approximations (|4.ip do not 
work at kd/n = 2, since this resonance is included in the 
neglected terms l/L[^ G ^ b ^ and 1/ L^' G ^^\ 

If R is larger than R cn , the square-gradient term in 
Eqs. (|2.2ip can be neglected. Then, after the substitu- 
tion of (|3.2ap into the first term, the mode-attenuation 



B. ASRB: Total and mode transmittances 

The preceding analysis of mode-attenuation lengths 
states that values of correlation length smaller than the 
first crossing point, R < R c x, result in both mode- 
attenuation lengths governed by intra-mode scattering 
only, see Eqs. (|4.ip . Therefore, mode and total trans- 
mittances are given by Eqs. (|3.8p and (|3.9p . however, 
with mode-backscattering lengths L]' and L 2 given by 
Eqs. (|2.21a[) and (|2.21cl) that are well described by their 
asymptotes (|4.1ap and (|4.1b[) . 

Also the preceding analysis says that values of cor- 
relation length larger than the second crossing point, 
R c2 < R, result in mode-attenuation lengths associated 
with inter-mode scattering only, see Eqs. (|4.2p ) and (|4.4p . 
In this case the theory that accounts for the transport 
properties should differ from that used in SRB waveg- 
uide. Within the approach presented here, the average 
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two-particle Green's function entering in Eq. (|2.3p . for 
the transmittance, is assumed to be equal to the product 
of two average one-particle ones. Therefore, to obtain the 
total transmittance (T(L)), we can substitute Eq. (|2.7[) 
into Eq. (|2 . 3f) . After evaluation of the integrals, the fol- 
lowing result becomes apparent 

N d 

(T(L)) = J2{T M (L/L n )), for R c2 < R. (4.5) 

n=l 

Here the total transmittance (T(L)) reads as a sum of av- 
erage mode transmittances (Tm{L / L n )) , describing the 
transparency of every nth propagating mode. Within the 
assumption necessary for the averaging procedure to be 
reasonable, they are determined by the following expres- 
sion 
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FIG. 5. (Color online) ASRB waveguide: Mode transmittance 
described by Eq. (|4.6p (solid line) and its asymptotes, given 
by (|4.7p and (|4.8|) for ballistic and diffusive regimes (dashed 
lines). The region where the diffusion terminates due to the 
influence of localization, is shown schematically. 



(T M (L/L n )) = 2 Ml-^ 



1 - exp ( - — 



(4.6) 



One can see that the mode transmittance is described 
by a function that depends on the parameter L/L n only. 
Note that in accordance with Eq. (|2.9p . the above pa- 
rameter depends on both backward and forward scatter- 
ing lengths explicitly given in Eqs. (|2.21j) and obeying 
the asymptotic expressions (|4.4[) . 

The mode transmittance (|4.6j) exhibits the ballistic be- 
havior for large mode attenuation length, 

(T M {L/L n ))nl-L/3L n for L < L n . (4.7) 

In this case the nth conducting channel is practically 
transparent. On the contrary, the mode transmittance 
is small when L n turns out to be much less than the 
length of the waveguide, 



(T M (L/L n ))n2L n /L for L n < L. 



(4.8) 



The inverse proportional dependence of Eq. (|4.8|) on the 
waveguide length L implies the diffusive regime of wave 
transport at a given nth mode. Figure [5] shows the be- 
havior of Eq. (|4~6| and its asymptotes (|4~7]) and (|4~8|) 
versus the parameter L/L n . 

It is necessary now to consolidate both aforementioned 
complementary approaches. Firstly, when R < R c i, the 
intra-mode scattering gives rise to the hierarchy of mode- 



backscattering lengths, < L { °' , see Eqs. (|4.1a[) and 
(|4.1bp . Due to this hierarchy and in accordance with 
two mode-transport regimes (|3.10j) and (|3.11j) . the to- 
tal transmittance (|3.9p can exhibit the ballistic, coex- 
istence and localized transport regimes. Secondly, when 
R C 2 < R, in accordance with two mode-transport regimes 
(|4.7p and (|4.8p , the total transmittance (|4.5[) can exhibit 
the ballistic and diffusive regimes. Here the coexistence 
transport regime can not emerge because the inter-mode 
scattering results in mode-attenuation lengths that ap- 
proximate each other, L\ « L2, see (|4.2b[) and (|4.4bl) . 
Specifically, the transport regimes of the total transmit- 
tance are defined as follows: 



(i) In the regime of localization, when the largest 

mode-backscattering length L b is smaller than the 
waveguide length L, and the correlation length R 
is smaller than the first crossing point R c i, 



1< L/L[ b) < L/L 2 ">, R<R cl , 



('■') 



(4.9) 



both propagating modes are strongly localized and 
their transmittances are exponentially small. The 
total transmittance is approximately equal to the 
transmittance of the first mode obeying the asymp- 
tote (|3TTT|) . 



(T(L)) « (T M (L/L[ b) )) 

7T 5 / 2 { L \ 



-3/2 



L 



16 



AL 



(6) 



exp 



AL 



(6) 



(4.10) 



with given by Eq. (|4.1ap . The waveguide is 
non-transparent in this regime. 

(ii) The coexistence regime arises when the smallest 
backscattering length i s smaller, while the 

largest backscattering length is larger than the 
waveguide length L, with the correlation length R 
smaller than the first crossing point R c i, 



L/L[ b) < 1< L/L ( 2 °>, R<R cl 



(b) 



(4.11) 



In this case, the first mode manifests the ballistic 
behavior p.lOp . while the second mode is still lo- 
calized in line with Eq. p. lip . Therefore, as before, 
the total transmittance (|3.9p is determined by the 
transmittance of the first mode, however governed 
by the ballistic asymptote (|3.10p . 



(T(L)) « (T M {L/Lf)) » 1 



(4.12) 



The coexistence of the ballistic and diffusive 
regimes is possible when the hierarchy of the mode- 



attenuation lengths, L2 < L^' , occurs 



-(6) 
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(iii) The diffusive regime emerges when both mode 
transmittances have the diffusive behavior and 
when the correlation length R is larger than the 
second crossing point R c2 , 



1 < L/Lx w L/L 2 , R c2 < R. 



(4.13) 



This regime arises if the attenuation is caused 
by the gradient-scattering mechanism. The mode 
transmittances have approximately the same be- 
havior that is described by Eq. (|4.8[) and, as a con- 
sequence, the total transmittance can be expressed 
as a double value of one of them, 



{T M {L/L{}) 
(T(L)) 



(T M (L/L 2 )), 
2{T M (L/L X )) 



(4.14a) 
4Lx/L. (4.14b) 



Here the approximate expression for L\ is given ei- 
ther by Eq. (@~2]) if kR c2 < kR < 1, or by Eq. (|4~4]) 
if 1 <C kR < fcA„. 

(iv) The ballistic regime emerges either under the con- 
ditions, 



L/L [b) < L/L 2 °> < 1 for R < R cl , 
L/Li w L/L 2 < 1 for R c2 < R, 



(4.15a) 
(4.15b) 

when the smallest mode-attenuation length, 
or L 2 , is larger than L. Both conducting chan- 
nels are open having approximately unit transmit- 
tances, (|3.10p or (|4.7[) . The total transmittance ap- 
proaches the total number of propagating modes, 



(T(L)) 



(4.16) 



Here the waveguide has almost perfect trans- 
parency. 

Note that our model is no able to describe the trans- 
mittances for a transition region that emerges within the 
interval R c \ < R < R c2 , when both intra- mode and inter- 
mode scattering have a comparable weight. However, 
this region is actually small in comparison with the large 
range of values of R for which the analytical approach 
has been developed. 



C. ASRB: Transmittance vs correlation length 

In order to discuss the dependence of the transmit- 
tances and their scaling parameters upon the dimension- 
less correlation length kR, Fig. ® is presented. In the 
region kR < kR c ±, where intra-mode scattering prevails, 
the transmittances (T M (L/ 'L,f>)) , (T M [L / L [ 2 b] ')) , (T(L)), 
the scaling parameters L/L\ \ L/L 2 b \ and the ratio of 

these parameters U b ' /L 2 \ are depicted. Otherwise, 
in the region kR c2 < kR, where inter-mode scattering 
prevails, the transmittances (Tm(L/Li)), (Tm{L/L 2 )), 
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FIG. 6. (Color online) ASRB waveguide: Mode- 
backscattering lengths, mode-attenuation lengths, Lh, , and 
transmittances, (Tj\/ n ) = {Tm (L / )) , (T(L)), vs dimen- 
sionless correlation length kR (here L„ = L„ for kR < kR c i 
and L„ = L n for kR C 2 < kR). The crossing points kR c i ~ 
0.13 and kR C 2 ~ 0.32 are indicated on the top of the upper 
plot. Note that there is no analytical results within the tran- 
sition region kR c \ < kR < kR C 2, in which the curves shown 
are obtained by simple interpolation. The global parameters 
are a/d = 0.0025 and kd/ir = 2.54, (ka = 0.02). In (a) the 
ratio L1/L2 is depicted; the asymptotes, L1/L2 ~ 35.78 and 
L1/L2 « 1 can be noted. In (b)-(c) the normalized waveguide 
length is L/d = 160; the boundary between the localized and 
coexistence regimes is at kR ~ 0.016, then, the boundary be- 
tween the coexistence and ballistic regimes is at kR ~ 0.05. 
In (d)-(e) L/d = 16, 000: the boundary between the localized 
and coexistence regimes is at kR ~ 0.075, then, the boundary 
between the diffusive and ballistic regimes is at kR « 7.0. 
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(T(L)), the scaling parameters L/Li, L/L2, and the ra- 
tio L1/L2, are depicted. In the first region the transmit- 
tances are computed from Eqs. (|3.9|) and (|3.8j) . whereas 
in the second region, from Eqs. (|4.5|) . (|4.6p . In both 
ranges, however, the scattering lengths are computed 
from Eqs. fOTjl . 

If i? decreases, the lengths L n and l[ 6 ' decreases as i? 3 
whereas the cycle length A„ has a fixed value, thus, the 
smallest value of R is restricted by the criterion (|2.15al) : 
here this criterion is fulfilled since Ai/2if ) w 0.03 < 1 

and A 2 /2L 2 b) rj 0.33 < 1 at the left boundary of the 
plots, kR = 0.01. Because the plots illustrate the case 
of weak correlations (R <C A„), the condition ()2.15b[) is 
satisfied automatically if (|2.15ap is met. 

The ratio of the lengths L± / j which is equal to the 
ratio L1/L2 for kR < kR c i, is in Fig. [5] (a). It exhibits 
the crossover from square-gradient scattering to gradient 
scattering. Therefore, with the increase of kR, the ratio 
decreases from the large value 



4 b) L X 



_±_ ~ -L ~ i 6 ( _i ) for kR < fci? cl , (4.17) 



to the unity, 
Li 



1 for fci? c2 < kR < fcA 2 . 



(4.18) 



The asymptotes ([4TT7|) and (|47T8]) arc obtained, respec- 
tively, through Eqs. (|4.1[) and (|4.4p . In contrast with the 
SRB waveguide, here the equality of mode-attenuation 
lengths, L\ w L 2 , emerges as a result of a strong in- 
tcrmode mixing, which is governed by the gradient- 
scattering mechanism. This non-trivial fact is in a strong 
contrast with a non-isotropic character of surface scatter- 
ing in the channel space [IH, H3] ■ 

In Fig.|5](b)-(c) and (d)-(e) two cases are illustrated in 
which, depending upon the roughness correlation length 
kR, different transport regimes can arise. The first case 
corresponds to a moderately large waveguide that can 
present the localized, coexistence or ballistic transport 
regime; see frames (b) and (c). In comparison, the second 
case corresponds to a larger waveguide that can present 
the diffusive transport regime in addition to the afore- 
mentioned ones; see frames (d) and (e). Note that differ- 
ent transport regimes, defined in Eqs. (|4.9[) . f|4.11[) . (|4.13j) 
and (|4.15|) . can be graphically identified with the use of 
frames (b) and (d) for the scaling parameters. 

The total and mode transmittances are plotted in 
Figs. [6] (c) and (e). These plots demonstrate that, due 

to the hierarchy of mode-backscattering lengths, < 



-(b) 



the transmittancc of the second mode is always 



smaller than the transmittancc of the first one, 



(T M {L/L 2 h) )) < (T M (L/L\°>)). 



(&)s 



(4.19) 



Therefore, within the localization (|4.9p and coexistence 
(14. lip regions, the total transmittance (|3.9p is mainly 



contributed by the first term, i.e. by the first-mode trans- 
mittance, 



(T(L))^(T M (L/L[ b) )). 



(4.20) 



Within the localization regime, (T(L)) behaves as indi- 
cated in Eq. (|4.10p . Within the coexistence regime the 
total transmittance is described by Eq. (|4.12p . In the 
diffusive regime the mode transmittances are equal, and 
behave as indicated in Eq. (|4.8p . The total transmittance 
has twice their value, see Fig.[5](e) and Eqs. (|4.14p . For 
the ballistic regime (|4.15bP , both conducting channels are 
open 



(Tm(L/L^)) 



(TM(L/L 2 b) )) 



1. 



(4.21) 



Here the approximation given by Eq. f)4.7[) describes well 
the behavior of the mode transmittances. Correspond- 
ingly, the value of the total transmittance tends to the 
total number N c i — 2 of the propagating modes, see 

Eq. ma . 

Note that the coexistence transport regime arises if 
the length L is such that L 2 & ' < L < , meanwhile the 



implied hierarchy < Lf' emerges if the correlation 
length is smaller than the crossing point, kR <C kR c \. 
Therefore, in the ASRB waveguide the scattering gov- 
erned by the square-gradient-scattcring mechanism is a 
necessary condition for the coexistence transport regime 
to occur. In comparison, if the gradient-scattering mech- 
anism prevails, the mode-attenuation lengths are almost 
equal, L\ rj L 2 . Then, with a continuous increment of 
L, the system passes directly from the ballistic to the 
diffusive regime. 



('<) 



D. ASRB: Transmittance vs mode parameter 

The mode-attenuation lengths and transmittances as 
functions of the mode parameter kd/ir are presented in 
Fig. [7] within the range of first two modes. The value of 
the parameters are selected in order to illustrate, within 
the range 2 < kdjix < 3, the coexistence transport regime 
of the ASRB waveguide, which is defined by Eq. (|4.1ip . 

Fig. [7] (a) shows the plots of L^/L and L { n' G)(h] / L. 
Here Ln^ is computed from Eqs. (|2.21a|) and (|2.21cl) . 
whereas L^ 4 ' -*^ is given by the first term of the above 
equations, i.e., 



1 



1 



-(A,G)(b) 



(A,G)(b) 



= 40.1 



(a/d^jhR) k 2 
(feid/7r) 2 (fc 2 d/7r) : 



■ exp 



-(fci 



! R 2 



If kd/ir increases, all the attenuation lengths also in- 
crease but exhibit characteristic dips at kd/ir = 1,2 
where a new normal mode opens. As we already noted, 
we can only indicate the position of these dips and do not 
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FIG. 7. (Color online) ASRB waveguide: Mode- 
backscattering lengths and transmittances versus the mode 
parameter for a/d = 0.0025, R/d « 0.005 and L/d = 160. 
Note that for fixed value of the correlation length R, and mode 
parameter within the interval 2 < kd/n < 3, the dimension- 
less correlation length kR ranges 0.031 kR ;$ 0.047, there- 
fore, the data correspond to the small-scale roughness (|3.ip 
and the coexistence regime (|4.11l) . 




describe their line shape in the narrow vicinity of integer 
values kd/n (see short discussion of these resonances in 
connection with Fig.@|. In contrast with the SRB waveg- 
uide, the curves depicting the first propagating mode of 
the ASRB waveguide show a dip at kd/n = 2. This 
dip is caused by the inter-mode scattering from the first 
propagating mode into the second, and this scattering 
is governed by the gradient-scattering mechanism. Re- 
markably, far from the resonances at kd/ir = 1,2, the 
lengths L{' G '^ b '/L and L 1 ^^^ / L are larger than the 

mode-backscattering lengths, / L and / L. This is 
because of the main contribution of the square-gradicnt- 
scattering lengths, Li /L and ° /L. There- 
fore, the approximate behavior of Ln/L is given by 
Eqs. (|4.1a|) and (|4.1b|) fnot plotted in the Figure). 



The total transmittancc p. 91) and mode transmit- 
tances (|3.8|) are plotted in Fig. 0(b). Note that the dips 
pertaining to the curves of L\ and Li become evident 
in the curves for mode transmittances (Tm(L/ L^)) and 

(Tm{L / L^)) ■ Certainly, the curve for the total trans- 
mittancc (T(L)) incorporates all the dips of the mode 
transmittances. Remarkably, without the weight of the 
square-gradient scattering mechanism, both mode trans- 
mittances would have the same value and the total trans- 
mittancc would have ballistic behavior. This fact can be 
graphically observed with the help of the curve related to 
the total transmittance given in Eq. (|3.22j) . in which the 
square-gradient-scattering mechanism has been ignored. 



V. CONCLUSIONS 

In this paper we report new results for transport prop- 
erties of quasi-lD surface-disordered waveguides (or elec- 
tronic conductors). Two models with symmetric and an- 
tisymmetric rough boundaries were studied in great de- 
tail. By examining these models, one passes from a sys- 
tem in which the attenuation of the propagating modes 
is associated with solely intra-mode scattering (the SRB 
waveguide), to a system in which both intra- and inter- 
mode scattering determine its attenuation (the ASRB 
waveguide). The comparison between the two models 
allows us to elucidate the role of all surface-scattering 
mechanisms. 

As was shown in our previous analytical studies [l7l — 
[2l| . one of mechanisms of scattering (the so-called 
square-gradient mechanism) is typically neglected in the 
literature since it emerges in the second order of pertur- 
bation theory. Specifically, the corresponding term in the 
expression for the inverse attenuation length is propor- 
tional to the squared variance of disorder, er 4 , and for 
this reason the square-gradient terms seem to be much 
smaller than the others describing conventional mecha- 
nisms of scattering (giving rise to the <r 2 -dependence). 
However, apart from the dependence on the disorder 
strength, one has to take into account the correlation 
properties of scattering profiles. It was found that the 
square-gradient mechanism essentially depends on the 
correlations of different nature as compared with the a 2 - 
terms in the total expression for the inverse attenuation 
length. 

In particular, it was argued that in some regions of pa- 
rameters, the square-gradient terms can prevail the a 2 - 
terms and, therefore, have to be considered as well, in 
spite of their apparent smallness. These predictions have 
been carefully analyzed in this paper, within two mod- 
els with symmetric and antisymmetric scattering profiles. 
Due to a principal difference between the two models, we 
analyze them separately. 

For the SRB waveguide the transmission occurs along 
two modes independently, due to the underlying symme- 
try of the model. Therefore, the transport through each 
of two modes can be described in terms of standard the- 
ories of localization in onc-dimcnsional models with ran- 
dom potentials. Thus, the total transmission is the sum 
of partial transmissions with the corresponding quantized 
transverse momentum. Our results demonstrate that in 
this case a typical situation occurs for which the trans- 
mission along two channels can be very different due to 
different values of corresponding backscattering lengths. 

For the ASRB waveguide there exist a competition of 
two surface-scattering mechanism, which in contrast with 
the SRB waveguide, induce the interplay between intra- 
mod and inter-mode scattering. This is a challenging sit- 
uation that persists with open questions. Nevertheless, 
the theoretical approach developed in Refs.[l8|,[2(| andl2ll 
allows one to determine intervals of parameters in which 
only one of the above mechanisms plays the main role, 
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i.e., it is possible to know when the channels can be con- 
sidered as two independent one-dimensional conductors 
and when they exhibit a strong interaction. Moreover, 
that theoretical approach indicates a rather small region 
of transition between the above possibilities, where ana- 
lytical results does not exist. The whole problem thus 
split off can be addressed by two complementary ap- 
proaches. When the intra- mode scattering prevails the 
transport through the waveguide can be described by 
the theory proposed for the SRB waveguide but with the 
proper scaling parameters. When the inter-mode scatter- 
ing prevails, it is necessary to apply the special approach 
developed here. Our results shows different transport 
regimes from localization to the ballistic, with the pos- 
sibility, depending upon the waveguide parameters, of 
finding the coexistence and diffusive regimes. The coex- 
istence regime can be present because of the hierarchy 



of mode-attenuation lengths associated with the square- 
gradient-scattering mechanism. Also, it is shown that 
the role of gradient mechanism of scattering is decisive, 
giving rise to a quite unexpected result. Specifically, it is 
demonstrated that in this case the partial transmissions 
along two channels are practically equal, and this effect is 
entirely due to the influence of the gradient terms in the 
expression for partial attenuation lengths. Due to the 
gradient-scattering mechanism that leads to the strong 
interaction between channels, the diffusive regime can 
take place. 

Our results can help to understand specific properties 
of surface scattering in waveguides and electronic devices, 
in dependence on correlation properties of scattering pro- 
files. In particular, one can understand the conditions 
under which the square-gradient mechanism of scatter- 
ing has to be taken into account. 
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